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Abstract 

A hierarchy of infinite-dimensional systems of hydrodynamic type 
is considered and a general scheme for classifying its reductions is 
provided. Wide families of integrable systems including, in partic- 
ular, those associated with energy-dependent spectral problems of 
Schrodinger type, are characterized as reductions of this hierarchy. 
iV-phase type reductions and their corresponding Dubrovin equations 
are analyzed. A symmetry transformation connecting different classes 
of reductions is formulated. 
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1 Introduction 



In this work we consider the system of evolution equations 
dY 

— =(A h Y), ^:=(A l r) + , *>0, (1) 

where (U, V) := UV X — U X V, and the function 

Y = Y{X,t), t:={to:=x,t u ...), 
is assumed to admit an expansion 

Y = H r 1 T75 1 , A — > OO. 

A A 

Here and henceforth we will denote by A = A(X) + + ^4(A)_ the standard 
decomposition of a power series of A in its components A + and A_ with pos- 
itive and strictly negative powers of A, respectively. The system of equations 
([]]) originated in the theory of the finite-gap solutions of the KdV equation 
(see [0) and, in a more general context, it appears in the analysis of the 
integrable hierarchies EDP(cf) (d > 1) associated with Schrodinger spectral 
problems with energy-dependent potentials M-M 

d-l 

d xx ^ = U(X, x)i/>, U(X, x) :=X d + Y, A*«i(x). (2) 

i=0 

The nonlinear evolution equations of the EDP(<i) hierarchy can be written 

as 

dtU = --A i>xxx + 2UA ijX + U x Ai, JU = (X l Y) + , (3) 
where the function Y can be determined by 

Y = 2Xl 

with ipi,ip2 being two independent solutions of (@). As a consequence of @ 
and it follows that 

2y xa; y - y x 2 - wy 2 + 4A d = o. (5) 
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This equation allows us to determine the coefficients of Y as differential 
polynomials depending on the potential coefficients {ui = Ui(x)}f~Q. In this 
way, the solutions of (|1|) provided by the EDP(d) hierarchy satisfy a reduction 
condition described by the differential constraint For example 

i) In the case of the EDP(l) hierarchy , the standard KdV hierarchy, we 

have U = A + uq and 

3 2 I 

u Q = -2y u y 2 = -y x + ^yi,**, ^ 
82/3 = - (yi,x) 2 + ^ViVi,xx + 82/12/2 + 2y 2 , xx . 

ii) The hierarchy EDP(2) is equivalent to the Zakharov-Shabat hierarchy || 

and one finds 

«i = -2yi, u = 3yl - 2y 2 , y 3 = V2y x y 2 - %y\ + y ljXX . (7) 

iii) For the case EDP(3) the first few relations arising from (gj) are 

u 2 = —22/1, Mi = 3yl - 2y 2 , u = -Ay\ + <oy x y 2 - 2y 3 . (8) 

The purpose of this paper is to study ([[]) as a hierarchy of nonlinear 
integrable models. In this sense, the form of (0) resembles that of the dis- 
persionless KP hierarchy [|Kj-[|nj] 

8Z 

— = {B h Z}, *>0, 

where {U, V} := U\V X — U X V\ is the standard Poisson bracket operation and 
Z = Z(X, t) is assumed to admit an expansion 

Zl (t) z 2 (t) 

Z = A H 1 — - H , A — ► 00. 

A A 

The subsequent analysis proves that the hierarchy (JI]), like the dispersion- 
less KP hierarchy, has deep connections with the theory of hydrodynamic 
systems. In fact, we notice that, in terms of the coefficients of the expan- 
sion of Y, the hierarchy (p]) becomes a set of infinite-dimensional systems of 
hydrodynamic type 

n 

diVn = ^2(y n -k,yi+k)i n > 1. (9) 
fe=l 
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We will show that the hierarchy exhibits a rich reduction theory which 
includes not only the standard types of reductions of the dispersionless KP 
hierarchy but also reductions of differential type which contain, in particular, 
the integrable hierarchies EDP(d). From this point of view the hierarchy flU) 
manifests a universal character. 

It should be observed that the zeros A = 7 of Y are Riemann invariants 
for (|l|) as they satisfy 

^ = Ml)d x j. (10) 

This fact allows one to write the dynamical equations of wide families of 
reductions in Riemann invariant form. In the case of differential reductions, 
the systems ([H]) have to be enlarged with appropriate equations of Dubrovin 

type0. 

One of the results of our work is the introduction of a symmetry transfor- 
mation for ([!]) which connects different classes of reductions. In particular, 
it transforms every EDP(d) hierarchy into its corresponding EDP(g? + 2) hi- 
erarchy. 



2 Reductions, TV-phase solutions and Dubrovin's 
equations 

2.1 General scheme 

Many interesting reductions of the hierarchy (|l|) can be formulated. We will 
be here concerned with those determined by constraints of the form 

F(\,Y o ,Y u ...,Y n )- = 0, Yj:=diY, (11) 

where it is assumed that 

dF 



^0. 



dY n 

Reductions of this type will be henceforth referred to as differential reductions 
of order n. The condition for (|TTJ) to determine one of such reductions is that 
for alH > 1 the function F must satisfy 



3=0 



j 
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provided Y verifies ([T]) and (PP. By taking into account that (|TI]) implies 

n 

(E F ^)- = ' 



j=0 

it follows that ([12]) is equivalent to 



n n 



(^/•^.'.l, - = o, (is) 



where 



a ' J - ./ J V j - 1 

Thus, it is obvious that a sufficient condition for ( |13"D to hold is that F 
satisfies a system of partial differential equations of the form 

n 

Fj-i - ^2 c kj F kY k +i-j = ajF + (3j, j = 1, . . . , n 

k=j ( 14 ) 

F n Y = a n+1 F + p n+1 , 
for a given set {ctj, /3j}?=i of entire functions of A (i.e. (aj)- = (/%)_ = 0). 

2.2 Differential reductions of order 

According to the general scheme, reductions of the form 

F(A,Y)_ = 0, (15) 

can be generated by solving 

dF 

Y w = aF + P, (16) 
for given functions a(A) and /5(A) satisfying «_=/?_ = 0. The general 



solution of ( 16) is 



a = 0, F = /5(A) InF + a (A), 

a *o, F = a( A)y*>-^, (17) 
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where a(A) is an arbitrary function of A. 

One of the simplest types of reductions ( |15D included in the class ([17]) is 

(\ N Y)_ = 0, N >1, (18) 

which means that Y is of the form 

r = i + ^ + ... + ^. (19) 

We will refer to the reduction ([19]) as the N -phase reduction of (|l|). Under the 
substitution (|I9| ) the hierarchy (|I]) becomes a finite system of hydrodynamic 
type for the coefficients {yi}f =1 . 

Example 

For 2-phase reductions 



r=1+ f + 5' (20) 



the system (|I|) becomes 



(21) 



dm = d x y 2 , diy 2 = (j/i,ya), 
d n yi = d n y 2 = 0, n > 2. 

Thus, by introducing a function Z such that 

yi = Z x , y 2 = Z t , t := tx, 

the 2-phase reductions of (|l|) are determined by the solutions of the homo- 
geneous Monge- Ampere equation 

Ztt = Z x Z tx — Z t Z xx , (22) 

which in turn reduces to the linear equation 

Wpp+pW pq + qW qq = 0, 

under the Legendre transformation (x, t) *—>■ (p, q) 

W = xp + tq - Z, p = Z x , q = Z t . 
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It should be noticed that an inhomogeneous version of the Monge-Ampere 
equation fl22|) describes the class of reductions of the dispersionless KP hier- 
archy which depend on two independent functions |T4|-||i6||. 

The iV-phase reductions admit a complete set of Riemann invariants sup- 
plied by the zeros {7;}^ of the function Y 

1 N 

i=l 

Thus, under the change of dependent variables 

y = (2/1, • • • , Vn) •-»• 7 = (7i, ■■■,1n), 
the system (|1|) can be written as 

Here (3i are the diagonal matrices 

{/3i) jk = M^ik = (7/ + 7f 1 2/i(7) + • • • + Vi{l))Sjk, 
the functions y n {~f) are given by the symmetric polynomials 

vnirt) = (-i) n X^i""^«> 

and the sum extends to all subsets J C {1, . . . , N} with n elements. Equa- 
tions (|23| ) are well known in the theory of iV-phase solutions of the KdV 
hierarchy |17| . 

2.3 Differential reductions of order 1 

Reductions of the form 

F(A,y J y JS )_ = o, (24) 

can be determined by solving 

dF dF 

Y W = aiF + (3l > Y w = 0i2F + ^ (25) 
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for given functions a i: Pi with (a»)_ = (A)- — 0- The compatibility condi- 
tions for these equations imply 

OC\ = — 1, «2 = 0, 

so that the general solution of the system (S3) is 



F = a(A)^ + ^l + 7 (A), (26) 

where a is an arbitrary function of A, and we are denoting a(X) := /^(A), 7(A) 
/?i(A). For a function F of the form (|26|) the reduction constraint ([24]) reads 

a(A)F x + a(A) =/5(A,t)F, (27) 

where 

Y x , a(A)' 



/3(A,t) := (a(A)^ + 



Y 

Examples 

1) The reduction corresponding to (p6|) with 

^=4 + p (28) 

where a is a given nonzero complex number, is determined by 

aY x + X = (X- yi )Y, (29) 

or, equivalently, by the following recurrence relation for the coefficients of 
the expansion of Y 

Vi+i = UiVi + ad x y { , i > 0. (30) 

2) If we set 

where a is a given nonzero complex number, then the reduction constraint 
pTj ) can be written as 

aY x + A 2 = (A 2 - \ yi + yi- y 2 )Y, (32) 
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which imposes the following recurrence relation 

Vi+2 = ym+i + (V2 - yl)Vi + atd x y h i > 0. (33) 

Under appropriate conditions it is possible to formulate combinations of 
both iV-phase and differential reductions. The coefficients {y n }^ =1 of the 
corresponding function Y are constrained by a certain system of differential 
equations with respect to the x variable. In terms of the zeros {^i}f = i of Y 
these differential equations are of Dubrovin type |TJ . 

For example, from (|2? ) it follows that the condition for differential reduc- 



tions of order 1 to admit a further iV-phase reduction is that the function 



N 



$:=\ N Y = H(\-~f t ), (34) 

satisfies 

a(A)$ z + \ N a{\) =/3(A,i)$. (35) 

This condition can be fulfilled only if (X N a)^ = 0. Furthermore, by inserting 
(|34] ) into (|35|) we get the following system of Dubrovin's equations 

** = >^— N - (36) 

Examples 



1) From ( p0|) we have that the 2-phase reduction of the differential reduction 
(|29|) is determined by 



2/2 = 2/i+ ad x y u y n = 0, n>3, (37) 

where y\ verifies the differential equation 

a 2 d xx y! + 3ayid x yi + yf = 0. (38) 

Let {7i}i =1 be the zeros of Y, then we can write yi = ji + 72, 2/2 = 7i72- 
Moreover, (|37D and (^) lead to the Dubrovin's equations 

ad xli = „ f r, i = 1,2. (39) 
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We observe that this reduction describes the 2-phase solutions of Burgers 
equation. Indeed, according to fl21|) and ( |37|) it follows 

u t = au xx + 2uu x , u := y u t := t x . (40) 

2) Equation (|33| ) implies that the 2-phase reduction satisfying the condition 
([32] ) is characterized by 

Vi = \{vl ~ «-^), Vn = } n>3, (41) 

where y\ verifies 

-2a 2 yi d xxyi + Aay\d xyi + 3a 2 (d xVl ) 2 - y\ = 0. (42) 
These equations lead to the Dubrovin's equations 

°^ = iufcri" ! ^ 2 ' (43) 

From fl21~D and (|4lf) we have that this reduction represents the 2-phase solu- 
tions of the evolution equation 

u t = --d x (a— ~ m2 )' u:=y!,t:= t v (44) 

2.4 Differential reductions of order 2 

Let us consider now differential reductions of second order 

F(\,Y,Y x ,Y xx )_ = 0. 

They are characterized by the solutions of 

dF dF 

Y w- Y -sY7r aiF+h ' 

dF dF 

Y w; +Y *wrr a2F+ ^ (45) 

dF 

Y— = a 3 F + f3 3 , 
oY x , 



■ XX 
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for given functions «j(A), /3j(A) with (a*)- = (/?»)- = 0. Equivalently, these 
conditions can be expressed as 

OF aiF+J3i Y 

^7 = p + yj("s^ + AO, 

dF a 2 F + ft Y „, R . .... 
<9F a 3 F + (3 3 



dY xx Y 

The compatibility conditions for these equations imply 

a t = -2, a 2 = a 3 = (3 2 = 0, 



and the following general solution of (46) arises 

^ = °(A)(^-^)+^+/*(A), (47) 

where a = a(A) is an arbitrary function and a; (A) := /3 3 (A), /3(A) := — /3 3 (A)/2. 
If we set 

a(A) = 2, a(A) = 4A d , (48) 

in (|47D, it follows that the corresponding reductions are the hierarchies EDP(d) 
associated to the Schrodinger spectral problems with energy- dependent po- 
tentials (H)-(0). Indeed it is enough to observe that fl4?D and (^) determine 
a reduction characterized by an equation of the form 

2Y XX Y - Yl - 4UY 2 + 4A d = 0, (49) 

where U is a polynomial in A 



U(\,t):=\ d + J2^(t), (50) 



i=0 



the coefficients of which can be recursively found from those of Y trough 
©■ 

From (pE9"D it is clear that the EDP(d) hierarchies admit iV-phase reduc- 
tions. If we introduce 



N 



$:=\ N Y = l[(\- ll ), (51) 



i=i 
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then (gU]) reads 

2$^$ - $2 - 4U$ 2 + AX d+2N = 0, (52) 
which is equivalent to the following system of Dubrovin's equations 

4 d+2N 

{dxlif = — T v»' i = l,...,N. (53) 

3 Symmetry transformations 

The hierarchy (|1|) exhibits an interesting type of symmetry transformations. 
Let 

Vi (*) 2/2 it) 
Y = 1 + ^ + ^ + ... 
A A 

be a solution of (P, and take a solution X = X(ti,t 2 , ■ ■ ■) of the associated 
system of differential equations 

d i X + y i (X,t 1 ,t 2 ,...) = 0, i>l, (54) 

then it follows that the function 

Y'(\,t'):=Y(\,t(t')), (55) 

where t' := (x', t[, t' 2 , . . . , t' n , . . . ) and 

t{t'):={X{x',t' l ,...),x',t' x ,... ,t' n -i^--)^ (56) 

is also a solution of (p. In order to prove this property we notice that from 
([]]) it follows at once that 

diAj - djAi = (A { , A,}, Ai := (X l Y) +J (57) 

which means that ([[]) defines a family of commuting vector fields in the space 
of Laurent series of A with coefficients depending on {yi}i>\. By setting A = 



in fl57j ) we find 

diyj - dffl = (yuyj). (58) 
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Hence, the contracted vector fields 

Vi := di - yid x , % > 1, (59) 

form a commutative family as well. In this way, given a solution Y = Y(X, t) 
of (0) then one has 

x? J r = ((A*- 1 r) + ,r> 1 , i>i, 

where (A, B)\ := Ad\B — Bd\A. Therefore, it readily follows that Y' is also a 
solution of (|l|) and, consequently, the transformation Y — > Y' is a symmetry 
of the hierarchy ([!]). 

Example 

An elementary solution of fl22|) is given by 

6 + 2 

It determines the following 2-phase solution of ([[]) 

y = 1 + T + — 

The associated system ( |54|) is 

<9X ,„ „ <9X 



— + bX = 0, ^ + ae" Wl 



which has the solution 



X = (c - at 2 )e Wl 



Hence, from fl55|) we get a new solution of ([I]) 

b(c — ati)e~ bx ae~ bx 
Y - 1+ A + — 

An interesting aspect of the transformation y — ► Y 7 is that it establishes 
some relationships between different differential reductions of (|l|). We notice 
that and (f>|) imply 

Y x = \(d x ,Y' + d x y 1 Y'), 

W , A (60) 

i^a = ^ (^'*' + d *yi Y x' + i^9 xx yi + d xx y 2 - y\d xx yi)Y J , 
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These formulas are useful to analyze the transformation properties of the 
differential reductions under Y — > Y' . For example, if Y satisfies a differential 
reduction (|27| ) of order 1 



a(X)Y x + a(X)=(3(X,t)Y, (61) 
then Y' satisfies a similar differential reduction 

a(X)Y^ + Xa(X) = (xp(X,t(t')) - a{X)d x y^Y'. (62) 

The way in which EDP(g?) reductions transform under fl5^|)-(]55|) is par- 
ticularly simple. Indeed, let us suppose that Y satisfies 

d-l 

2Y XX Y - Y 2 - AUY 2 + AX d = 0, U(X, x) := X d + A ^(z), ( 63 ) 

i=0 

then, as a consequence of the transformed function Y' verifies 

1Y**Y' ~ {Y^f ~ WY' 2 + 4A d+2 = 0, (64) 

where 

U' := X 2 U - -(Xd^yx + d xx y 2 - yid xx yx) + ^{d x yi) 2 . (65) 

This means that Y — > Y' transforms the EDP(cf) hierarchy into the EDP(g? + 
2) hierarchy. In particular, this result proves that the whole family of hier- 
archies of integrable models associated with energy-dependent Schrodinger 
problems can be generated from its two first members, namely: the KdV 
hierarchy and the Zakharov-Shabat hierarchy. 

From (^) one can also derive the transformation properties of Dubrovin's 
equations under the symmetry Y — > Y' . Obviously, if Y satisfies a iV-phase 
reduction so does Y'. Moreover, the zeros of Y transform as 7^(t') = ji(t(t')) 
and (|60D lead to identities of the type d^j^ = ^pxli- 
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